------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  ...\1903 JOP\trust_poa.log
  log type:  text
 opened on:  12 Mar 2019, 12:40:00

. 
. set seed X6ebfb33df1c34e2040655e83559d6291000412ea

. 
. logit trust_poa pattern reform ///
>      saccounty sac_pattern sac_reform ///
>          black black_pattern black_reform

Iteration 0:   log likelihood = -1603.8814  
Iteration 1:   log likelihood = -1583.1375  
Iteration 2:   log likelihood = -1583.0527  
Iteration 3:   log likelihood = -1583.0527  

Logistic regression                               Number of obs   =       2323
                                                  LR chi2(8)      =      41.66
                                                  Prob > chi2     =     0.0000
Log likelihood = -1583.0527                       Pseudo R2       =     0.0130

-------------------------------------------------------------------------------
    trust_poa |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
      pattern |   .0302732   .1455286     0.21   0.835    -.2549577    .3155041
       reform |   .1614548   .1429628     1.13   0.259    -.1187471    .4416568
    saccounty |   .0536775    .141655     0.38   0.705    -.2239612    .3313163
  sac_pattern |   -.046725   .2043065    -0.23   0.819    -.4471583    .3537083
   sac_reform |   .0176118   .2053752     0.09   0.932    -.3849163    .4201399
        black |   -.705356   .2604185    -2.71   0.007    -1.215767   -.1949451
black_pattern |   .0144719   .3823897     0.04   0.970    -.7349982    .7639419
 black_reform |  -.6722452   .3974059    -1.69   0.091    -1.451146     .106656
        _cons |  -.1528885   .0998926    -1.53   0.126    -.3486744    .0428973
-------------------------------------------------------------------------------

. 
. estsimp logit trust_poa pattern reform ///
>      saccounty sac_pattern sac_reform ///
>          black black_pattern black_reform

Iteration 0:   log likelihood = -1603.8814
Iteration 1:   log likelihood = -1583.2394
Iteration 2:   log likelihood =  -1583.053
Iteration 3:   log likelihood = -1583.0527

Logistic regression                               Number of obs   =       2323
                                                  LR chi2(8)      =      41.66
                                                  Prob > chi2     =     0.0000
Log likelihood = -1583.0527                       Pseudo R2       =     0.0130

------------------------------------------------------------------------------
   trust_poa |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     pattern |   .0302732   .1455286     0.21   0.835    -.2549577    .3155041
      reform |   .1614548   .1429628     1.13   0.259    -.1187471    .4416568
   saccounty |   .0536775    .141655     0.38   0.705    -.2239612    .3313163
 sac_pattern |   -.046725   .2043065    -0.23   0.819    -.4471583    .3537083
  sac_reform |   .0176118   .2053752     0.09   0.932    -.3849162    .4201399
       black |   -.705356   .2604185    -2.71   0.007    -1.215767   -.1949451
black_patt~n |   .0144719   .3823897     0.04   0.970    -.7349982    .7639419
black_reform |  -.6722452    .397405    -1.69   0.091    -1.451145    .1066543
       _cons |  -.1528885   .0998926    -1.53   0.126    -.3486744    .0428973
------------------------------------------------------------------------------

Simulating main parameters.  Please wait....
% of simulations completed: 11% 22% 33% 44% 55% 66% 77% 88% 100% 

Number of simulations  : 1000
Names of new variables : b1 b2 b3 b4 b5 b6 b7 b8 b9

. 
. 
. * Control group, non-Sac
. 
. setx median

. simqi, prval(1) genpr(control_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         0       median  
       reform |         0       median  
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~a=1) |   .4616325      .024781     .4163793    .5105217

Simqi generated the following new variable(s): control_pr

. 
. simqi, fd(prval(1) genpr(pattern_fd)) changex(pattern 0 1)

First Difference: pattern 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~a = 1) |   .0073649     .0359705    -.0648496     .081571

Simqi generated the following new variable(s): pattern_fd

. simqi, fd(prval(1) genpr(reform_fd)) changex(reform 0 1)

First Difference: reform 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~a = 1) |   .0396602     .0347875    -.0283506    .1066135

Simqi generated the following new variable(s): reform_fd

. 
. 
. simqi, fd(prval(1) genpr(sac_fd)) changex(saccounty 0 1)

First Difference: saccounty 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~a = 1) |   .0129925     .0359387    -.0578316    .0862136

Simqi generated the following new variable(s): sac_fd

. simqi, fd(prval(1) genpr(black_fd)) changex(black 0 1)

First Difference: black 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~a = 1) |  -.1618542      .056247    -.2697556   -.0491709

Simqi generated the following new variable(s): black_fd

. 
. 
. * Pattern group, non-Sac
. 
. setx pattern 1

. simqi, prval(1) genpr(pattern_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         1         1     
       reform |         0       median  
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~a=1) |   .4689974     .0261953     .4191036     .523754

Simqi generated the following new variable(s): pattern_pr

. 
. 
. * Reform group, non-Sac
. 
. setx pattern 0

. setx reform 1

. simqi, prval(1) genpr(reform_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         0         0     
       reform |         1         1     
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~a=1) |   .5012926     .0246626     .4551874    .5520061

Simqi generated the following new variable(s): reform_pr

. 
. 
. * Difference in baseline probabilities
. 
. tabstat control_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
  control_pr |  .4616325   .024781  .3737745  .4213502  .4286569  .4621551  .4925119  .5012761  .5292172
--------------------------------------------------------------------------------------------------------

. tabstat pattern_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
  pattern_pr |  .4689974  .0261953   .381394  .4260152  .4365687  .4676998  .5028521  .5165344  .5561548
--------------------------------------------------------------------------------------------------------

. tabstat reform_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
   reform_pr |  .5012926  .0246626   .417805   .462192  .4706634  .5000845  .5340519  .5430182  .5738892
--------------------------------------------------------------------------------------------------------

. 
. gen dif_con_pat_pr = pattern_pr - control_pr
(1759 missing values generated)

. gen dif_con_ref_pr = reform_pr - control_pr
(1759 missing values generated)

. gen dif_pat_ref_pr = reform_pr - pattern_pr
(1759 missing values generated)

. 
. tabstat dif_con_pat_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_con_pa~r |  .0073649  .0359705 -.0915518 -.0504381 -.0393936  .0071583   .054453  .0680511  .1155621
--------------------------------------------------------------------------------------------------------

. tabstat dif_con_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_con_re~r |  .0396602  .0347875 -.0683937 -.0174947 -.0062365  .0406145  .0850877  .0977725   .168308
--------------------------------------------------------------------------------------------------------

. tabstat dif_pat_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_pat_re~r |  .0322953  .0357313  -.099369 -.0252049  -.011191  .0328705  .0770098  .0916809  .1571268
--------------------------------------------------------------------------------------------------------

. 
. 
. * Control group, Sac
. 
. setx median

. setx saccounty 1

. simqi, prval(1) genpr(sac_control_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         0       median  
       reform |         0       median  
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         1         1     
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~a=1) |    .474625     .0266797     .4213253    .5278913

Simqi generated the following new variable(s): sac_control_pr

. 
. simqi, fd(prval(1) genpr(sac_pattern_fd)) changex(pattern 0 1 sac_pattern 0 1)

First Difference: pattern 0 1 sac_pattern 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~a = 1) |  -.0016983     .0367534    -.0734718    .0718328

Simqi generated the following new variable(s): sac_pattern_fd

. simqi, fd(prval(1) genpr(sac_reform_fd)) changex(reform 0 1 sac_reform 0 1)

First Difference: reform 0 1 sac_reform 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~a = 1) |   .0442381     .0381929    -.0315103    .1224265

Simqi generated the following new variable(s): sac_reform_fd

. 
. 
. * Pattern group, Sac
. 
. setx pattern 1

. setx sac_pattern 1

. simqi, prval(1) genpr(sac_pattern_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         1         1     
       reform |         0       median  
  sac_pattern |         1         1     
   sac_reform |         0       median  
    saccounty |         1         1     
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~a=1) |   .4729266     .0262612     .4224745    .5261227

Simqi generated the following new variable(s): sac_pattern_pr

. 
. 
. * Reform group, Sac
. 
. setx pattern 0

. setx sac_pattern 0

. setx reform 1

. setx sac_reform 1

. simqi, prval(1) genpr(sac_reform_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         0         0     
       reform |         1         1     
  sac_pattern |         0         0     
   sac_reform |         1         1     
    saccounty |         1         1     
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~a=1) |   .5188631     .0270311     .4646192    .5703435

Simqi generated the following new variable(s): sac_reform_pr

. 
. 
. * Difference in baseline probabilities
. 
. tabstat sac_control_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
sac_contro~r |   .474625  .0266797  .3953791  .4312263  .4411028  .4741347  .5104108  .5197653  .5532965
--------------------------------------------------------------------------------------------------------

. tabstat sac_pattern_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
sac_patter~r |  .4729266  .0262612  .4030392  .4307849  .4408427   .472029  .5069598  .5182082  .5648054
--------------------------------------------------------------------------------------------------------

. tabstat sac_reform_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
sac_reform~r |  .5188631  .0270311  .4392072  .4731504  .4843222  .5190101  .5543062  .5618708  .6073281
--------------------------------------------------------------------------------------------------------

. 
. gen dif_sac_con_pat_pr = sac_pattern_pr - sac_control_pr
(1759 missing values generated)

. gen dif_sac_con_ref_pr = sac_reform_pr - sac_control_pr
(1759 missing values generated)

. gen dif_sac_pat_ref_pr = sac_reform_pr - sac_pattern_pr
(1759 missing values generated)

. 
. tabstat dif_sac_con_pat_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_sac~t_pr | -.0016983  .0367534 -.1153741 -.0638004 -.0472047 -.0019202  .0457835  .0576247  .1139996
--------------------------------------------------------------------------------------------------------

. tabstat dif_sac_con_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_sac_co.. |  .0442381  .0381929 -.0871337 -.0179733 -.0044478  .0435036  .0944596  .1070256  .1758039
--------------------------------------------------------------------------------------------------------

. tabstat dif_sac_pat_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_sac_pa~r |  .0459364  .0362748 -.0660989 -.0141295 -.0024831  .0480253  .0906137  .1024147  .1520375
--------------------------------------------------------------------------------------------------------

. 
. 
. * Control group, Black
. 
. setx median

. setx black 1

. simqi, prval(1) genpr(black_control_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         1         1     
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         0       median  
       reform |         0       median  
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~a=1) |   .2997782     .0537236     .2056278    .4092354

Simqi generated the following new variable(s): black_control_pr

. 
. simqi, fd(prval(1) genpr(black_pattern_fd)) changex(pattern 0 1 black_pattern 0 1)

First Difference: pattern 0 1 black_pattern 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~a = 1) |   .0094619     .0790728    -.1441042    .1666322

Simqi generated the following new variable(s): black_pattern_fd

. simqi, fd(prval(1) genpr(black_reform_fd)) changex(reform 0 1 black_reform 0 1)

First Difference: reform 0 1 black_reform 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~a = 1) |  -.0918135     .0704394    -.2280496    .0438154

Simqi generated the following new variable(s): black_reform_fd

. 
. 
. * Pattern group, Black
. 
. setx pattern 1

. setx black_pattern 1

. simqi, prval(1) genpr(black_pattern_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         1         1     
black_pattern |         1         1     
 black_reform |         0       median  
      pattern |         1         1     
       reform |         0       median  
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~a=1) |   .3092401     .0581914     .2073215    .4247012

Simqi generated the following new variable(s): black_pattern_pr

. 
. 
. * Reform group, Black
. 
. setx pattern 0

. setx black_pattern 0

. setx reform 1

. setx black_reform 1

. simqi, prval(1) genpr(black_reform_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         1         1     
black_pattern |         0         0     
 black_reform |         1         1     
      pattern |         0         0     
       reform |         1         1     
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~a=1) |   .2079647     .0471983      .122954    .3044373

Simqi generated the following new variable(s): black_reform_pr

. 
. 
. * Difference in baseline probabilities
. 
. tabstat black_control_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
black_cont~r |  .2997782  .0537236  .1579724  .2176589  .2323455  .2957931  .3707047  .3943076  .4667004
--------------------------------------------------------------------------------------------------------

. tabstat black_pattern_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
black_patt~r |  .3092401  .0581914  .1488188  .2180398  .2348376  .3062351  .3836353  .4058624  .5772182
--------------------------------------------------------------------------------------------------------

. tabstat black_reform_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
black_refo~r |  .2079647  .0471983  .0961624  .1327485  .1491737  .2061257  .2697031  .2867158  .4037836
--------------------------------------------------------------------------------------------------------

. 
. gen dif_black_con_pat_pr = black_pattern_pr - black_control_pr
(1759 missing values generated)

. gen dif_black_con_ref_pr = black_reform_pr - black_control_pr
(1759 missing values generated)

. gen dif_black_pat_ref_pr = black_reform_pr - black_pattern_pr
(1759 missing values generated)

. 
. tabstat dif_black_con_pat_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_bla~t_pr |  .0094619  .0790728 -.2332804 -.1247284 -.0909142  .0106966   .109612  .1411162  .2907198
--------------------------------------------------------------------------------------------------------

. tabstat dif_black_con_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_black_.. | -.0918135  .0704394 -.3019089 -.2118179 -.1863694 -.0910894 -.0064056   .022456  .1144115
--------------------------------------------------------------------------------------------------------

. tabstat dif_black_pat_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_black_p~ | -.1012753  .0738581 -.4526429 -.2250232 -.1934426 -.1003242 -.0092614   .024904  .1622899
--------------------------------------------------------------------------------------------------------

. 
. 
. drop b1-b9

. drop control_pr-dif_black_pat_ref_pr

. 
. log close
      name:  <unnamed>
       log:  ...\1903 JOP\trust_poa.log
  log type:  text
 closed on:  12 Mar 2019, 12:40:01
------------------------------------------------------------------------------------------------------------
